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ABSTRACT 

A fraction of the first generation of stars in the early Universe may be very massive (> 300 Mq) 
as they form in metal-free environments. Formation of black holes from these stars can be 
accompanied by supermassive collapsars to produce long gamma-ray bursts of a unique type 
having a very high total energy 10®"^ erg) as recently suggested by several authors. We 
present new stellar evolution models of very massive Population III stars including the effect 
of rotation to provide theoretical constraints on super-collapsar progenitors. We find that the 
angular momentum condition for super-collapsar can be fulfilled if magnetic torques are ignored, 
in which case Eddington-Sweet circulations play the dominant role for the transport of angular 
momentum. We further find that the initial mass range for super-collapsar progenitors would be 
limited to 300 Mq ^ M < 700 Mq. However, all of our very massive star models of this mass 
range end their lives as red supergiants rather than blue supergiants, in good agreement with 
most of the previous studies. The predicted final fate of these stars is either a jet-powered type 
IIP supernova or an ultra-long, relatively faint gamma-ray transient, depending on the initial 
amount of angular momentum. 


Subject headings: stars: evolution — stars: rotation — stars: Population III — gamma-ray bursts — 
early universe 


1. Introduction 


Population III stars are the main sources of ion¬ 
izing photons and the only producers of heavy el¬ 
ements in the early Universe. Their exact role 
on the reionization history and the chemical evo¬ 
lution largely depends on their initial mass func¬ 
tions. In the metal-free environment where the 
first stars form, thermal cooling of star-forming re¬ 
gions is dominated by hydrogen molecules, which 
are much less efficient coolants than heavy ele¬ 
ments and dusts. It is therefore widely believed 
that Population III (Pop III) stars are system¬ 
atically m ore massive than Pop II and Pop I 
stars (e.g.. lAbel et al.l[2002l : iBroinm et al. 2002 ). 

According to recent numerical simulations, 
most of Pop III stars would have a mass range 
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of 10 - 1000 M 

P) (e 

g.. Omukai & Palla 20031 

Ohkubo et al.l2009l: Bromm et al.ll2009l: Hosokawa et al. 

2011: Stacv et al. 

2012 

: Hirano et al. 2014 

), while 


formation of more massive Pop III stars (> 


ITOO Mq) still remains a good possibility (IHosokawa et al 


20121) . A significant fraction of them would pro¬ 


duce core-collapse supernovae for initial masses 
of 10 — 100 Mq, and pulsational pair-instability 
or pair-instability supernovae for 100 — 300 Mf^ 

( Heger fc Wooslet^l2002l : Umeda fc Nomotoll2002l) . 

According to the recent numerical simulations, the 
first stars may ro tate at a velocity close to the 
critic al rotation ( Stacv et al. 2011 : Stay et al~l 
l2013l) . and these mass ranges can be signifi¬ 
cantly lowered if rapid rotati on induces strong 
chemical mixing inside stars (lYoon et al.l 1201 


ChatzoDOulos fc Wheeler! 2012 : cf. Glatzel et al 

19851 ). Gamma-rav bursts may also occur with 
sufficiently high ini tial angular mom entum, for 


about 12 — 80 Mq (lYoon et al.ll2012l) . Detection 
of these events at z ~ 20 with next generations 
of telescopes will give us invaluable information 
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on the nature of Po 

p III stars fe.g.. Kawai et al. 

2006: Greiner et al. 

2009: Salvaterra et al.l 20091: 

Tanvir et al. 

2 OO 9 I: ICucchiara et al.l 201lll. 


On the other hand, Pop III stars with initial 
masses of 300 — 50000 Mq are supp osed to di¬ 
rectly collapse to a black hole (e.g., iFrver et al 


2 OOII : iHeeer fc WooslevI 1200211 . In particular, 


those with Mjnit > 1000 Mq are often consid¬ 
ered seeds for super-massiv e black holes that are 
found at high redshift (e.g. iMadau fc Reed 2001 


ISchneider et al. 2002ll . Recently several authors 
suggested that formation of black holes in very 
massive stars can be accompanied by a super- 
collapsar (i.e, collapsar with a black hole mass 
higher than a few hundr eds solar masses) to 
prod u ce relativistic jets (iKomissarov fc Barkov 


2010 ; 


Meszaros fc Rees 201ot ^wa fc lokal 201ll; 


Maio fc Barkov 2014ll . Given that the considered 


black holes are very massive compared to the case 
of ordinary gamma-ray bursts, there might be 
unique observational signatures from these super- 
collapars, depending on the details of the progen¬ 
itor structure. For example, the amount of the 
total energy released by jets could be significantly 
l arger than those of ordinary gamma-ray bursts 
(Komissarov fc Barkov! l2010t Meszaros fc ReesI 


201 ( 111 . 


One of the key conditions for having a collap¬ 
sar event is rapid rotation of the infalling matter 
onto the black hole, of which the specific angu¬ 
lar momentum must be high enough to form an 
accretion disk (i.e., j > 10^® — 10^® c m^ s for 
the black hole mass of 3 to 1000 Mq; IWooslev 


1993 : MacFadven fc Wooslev 1999ll : If the black 


hole mass is higher than a few hundreds solar 
masses, the corresponding radius of the disk is too 
large for ne utrino annihilation to produc e a rela¬ 
tivistic jet ( MacFadven fc Wooslev 1999l l. There¬ 
fore, another key ingredient for super-collapsar is 
strong magnetic fields (i.e., B > 10^^ G around 
the black hole) such that the rotational energy of 
the black hole may be extracted by the Blanford- 
Znajek mechanism (IKomissarov fc Barkov 201(1 


Meszaros fc ReeJ201(ll l. The progenitors of super- 

collapsars should be both rapidly rotating and 
strongly magnetic. 

IYoou et al.l ( 20121 hearafter, YDL12) calcu¬ 
lated evolutionary models of Pop III stars includ¬ 
ing rotation and magnetic fields. In their mod¬ 


els the Tayler-Spruit (TS) dynamo (jSDruitJ 1200211 


plays a key role in transporting angular momen¬ 
tum inside stars. Because of the efficient angular 
momentum transport from the convective core to 
the outer layers via magnetic torques, stars can 
maintain near-rigid rotation throughout the main 
sequence, and lose angular momentum via mass 
shedding as the surface reaches critical rotation. 
This effect becomes more important for a very 
massive Pop III star of which the surface luminos¬ 
ity is close to the Eddington limit because it can 
reach the critical rotation even with a fairly low 
rotation velocity. As a result. Pop III stars more 
massive than about 300 Mq cannot retain enough 
angular momentum to produce a super-collapsar, 
according to YDL12. 


The TS dynamo adopted by YDL12 is one of 
the most considered explanations for the rapid 
transport of angular momentum in radiative lay¬ 
ers of stars, which is implied by a number of 
obser vations, in particu l ar wi t h low-mass stars 


e.g.. lEggenbereer et al 


20051: ISniis et al.l 120081: 


Eggenberger et al. 2012 : Cantiello et al.l 201 41 . 


However, the efficiency of the TS dyn amo is still a 
matt e r of debate in the community ( Braithwaitel 
2006 : ^hn et al.l 200?! : IPenissenkoy et alT 2010t 
Arlt fc Riidigeij 201 lb . The strong braking im¬ 


plied in low-mass stars might be instead related 
to so me other mechanisms like internal gravity 
waves ( Zahn et al. 1997 ; I Talon et al. 2002 ), which 
have not been well understood in the context of 
massive stars and thus might not be necessarily 
relevant. If the transport of angular momentum 
in massive stars is much less efficient than what is 
predicted with the TS dynamo, they could retain 
a large amount of angular momentum inside the 
core, while magnetic fields of large scales might 
be generated in convective layers before collapse 
or in the accretion disk after black hole formation, 
thus fulfilling the necessary conditions for super- 
collapsar. Given the importance of super-collapsar 
events as a probe on the nature of Pop III stars, it 
is worthwhile to investigate this possibility in some 
detail to provide a theoretical boundary condition 
for super-collapsar progenitors. For this purpose, 
we present new evolutionary models of very mas¬ 
sive Pop III stars (Minit ^ 300 Mq) without mag¬ 
netic torques according to the TS dynamo. This 
means that we only consider Eddington-Sweet cir¬ 
culations and other hydrodynamic instabilities for 
the transport of angular momentum. 
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We briefly explain the numerical methods used 
for our models in Sect. [2l In the section that 
follows (Sect. |3]), we describe the properties of 
very massive Pop III stars on the zero-age main 
sequence (ZAMS). In particular, we show that 
the angular momentum condition for collapsar be¬ 
comes more difficult to meet for a higher initial 
mass because the Eddington factor at the stellar 
surface becomes larger. In Sect. IH we present the 
evolutionary models including rotation and dis¬ 
cuss the initial conditions of Pop III stars needed 
for super-collapsar progenitors. We discuss the im¬ 
plications of our results for the final fate of super- 
collapsar progenitors and the relevant uncertain¬ 
ties of our models in Sect. [S] We conclude our 
discussions in Sect. El 


2. Numerical methods 


We used the one-dimensional hydrodynamic 
stell ar evolution cod e that is described in YDL12 
and iKozvreva et al. ( 2014 ). This code implic¬ 
itly solves the stellar structure equations for 
which the Henyey-typ e method is adopted. T he 
OPAL opacity table ( Isrlesias fc Rogers 1990ll is 
used for T > 10"^ K and the prescription by 
Alexander fc Fer gusonl (I l994l) for low er tempera¬ 
tures. We follow lEndal fc Sofial (|1976ll to consider 
the effect of the centrifugal for ce o n the stellar 
scructure, and Heger et al. ( 200nh and iHeger et al 


l)2005[ ) for the redistribution of angular momen¬ 
tum, which is approximated as a diffusive process. 
We used the same physical parameters for con¬ 
vection, semi-convection and rotationally-induced 
chemical mixing and transport of angular momen¬ 
tum, but the TS dynamo is not included in the 
present study except for one model sequence. This 
means that the considered transport processes due 
to rotation include Eddington-Sweet circulations, 
the secular shear instability and the Goldreich, 
S chubert and Ericke (GSE) instability as explained 
in iHeger et al. ( 2000l ). The overshooting is applied 
for 0.335 times local pressure scale heights above 
the convectively unstable core and the adopted 


semi- convection para meter (assmi ; Langer et al 


1983h is 1.0, following iBrott et al. ( 201lh . Com¬ 
pared to YDL12, the nuclear network was im- 
proved to trea t silico n bu rning as d i scusse d in 
Kozyreva et al. ( 2014 ) and Kozyrev^ ( 2014 ). In 
some of our models, the silicon burning was fol¬ 
lowed with 13 main isotopes. In this way the main 


alpha-chain can be described without considering 
neutronization that does not play an important 
role in the structure of very massive stars under¬ 
going the pair-instability. 

The treatment of mass loss for Pop III stars 
and its physical justification are fully described in 
YDL12. In short, the mass loss in our models is 
dominated by the centrifugally-driven mass shed¬ 
ding, which is treated as the following: 


M(z;)=max 10 M 0 yr ^,M(u = 0) 
where 


1 


1 - n 
( 1 ) 


0.43 


:=- and Ucrit = 


GM 


(l-T). (2) 


The critical rotation r^crit represents the modi¬ 
fied Keplerian limit considering the effect of ra¬ 
diation when the stars approach the Edding¬ 
ton limit, for which the Eddington factor T is 
given by kL/{4:'kGMc). M(v = 0) is the wind 
mass loss rate for the non-rotating case , for 
which we followed Kudritzki et al.l (119891) and 
Nieuwenhuiizen fc de JageiT i 1990l) for > 10"^ 
and Teff < 10"^, respectively, with metallicity de¬ 
pendence of Here, the metallicity means 

the total mass fraction of CNO elements at the 
surface, and in this way the enhancement of the 
mass loss rate resulting from surface enrichment 
of heavy elements by chemical mixing is taken into 
account (See YDL12 for more discussion on this 
issue). 


3. Properties on the zero-age main-sequence 

Before presenting the evolutionary models, here 
we discuss the properties of very massive Pop III 
stars on the ZAMS. In Table 1, we provide the in¬ 
formation on the ZAMS models for the mass range 
of 300 - 20000 Mq. Here we assume that stars on 
the ZAMS rotate as a solid body, which can be 
justified because of the expected rapid transport 
of angular momentum in chemically homogeneous 
stars by con vection and Eddington -Sweet circu¬ 
lations fe.g.. iHaemmerle et al.l [20131 see also the 
discussion below). 

The surface luminosities of these very massive 
stars are close to the Eddington limit: the Edding¬ 
ton factor varies from P = 0.774 at 300 Mq to P = 
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0.996 at 20000 Mq. In consequence, these stars 
cannot have a very high critical rotation veloc¬ 
ity, and the total amounts of angular momentum 
should be limited accordingly. The corresponding 
^'crit/'^K (= vi - T; see Eq. with solid-body ro¬ 
tation changes from 0.48 to 0.06. Because the stel¬ 
lar radius is larger with higher mass, both the total 
angular momentum and the specific angular mo¬ 
mentum at the critical rotation gradually increase 
with increasing mass up to about 10000 M 0 , de¬ 
spite the fact that Vait/vK decreases. However, 
VcTit/vK becomes so small for 20000 Mq that the 
specific angular momentum becomes smaller than 
that of 10000 Mq star. Another factor that limits 
the amount of angular momentum is the change of 
stellar structure as T approaches zero. The dimen¬ 
sionless radius of gyration (i.e., k = J/ {vrotMR); 
see Table 1) decreases with increasing mass and 
therefore the amount of angular momentum for a 
given set of Ucrit, M and R should decrease as well. 


This property leads to the following important 
conclusion: Pop III stars with Mi > 3000 Mq 
at their birth cannot contain enough angular mo¬ 
mentum in their cores to produce a collapsar be¬ 
cause of the effect of high radiation pressure, as 
illustrated in Fig. [T] Only the outermost layers 
have specific angular momentum above the crit¬ 
ical limit in these stars, and therefore produc¬ 
tion of an energetic gamma-ray burst cannot occur 
(see discussion below). This fact is based only on 
the structure of Pop III stars on the ZAMS, and 
does not depend on the uncertain physical pro¬ 
cesses like mass loss that has strong impact on the 
stellar evolution. This gives a strong upper mass 
limit for super-collapsar progenitors. Given that 
Mi > 300 Mq is required to avoid a pair-instability 
explosion, only Pop III stars in the mass range of 
300 Mq < Mi < 3000 Mq can be considered po¬ 
tential super-collapsar progenitors that can pro¬ 
duce an energetic gamma-ray burst, in terms of 
the initial angular momentum budget. Therefore, 
this result effectively rules out the possibility of 
an energetic gamma-ray burst with formation of 
intermediate-mass black holes of ^ 10"^ Mq from 
super-massive Pop III stars that can be seeds for 
sup er-massive black ho l es in high-redshift q uasars 
(cf. ISoolvar et al. 20091 : Freese et al. 2010ll . 


4. Evolutionary models 

We calculated evolutionary models for 300, 500, 
700, 1000, and 2000 Mq including rotation but 
without the TS dynamo as summarized in Ta¬ 
ble. 2. We exceptively included the TS dynamo 
in the sequence S300M that shares the same ini¬ 
tial condition with S300C, for comparison. It is 
well known that these very massive stars undergo 
the pair instability when the central temperature 
exceeds about 10® K, which leads to rapid con¬ 
traction of the core and the consequent rapid oxy¬ 
gen and silicon burning on a dynamical timescale. 
Because of the very high binding energy, this ex¬ 
plosive nuclear burning cannot make these stars 
explode to produce a supernova. This pair insta¬ 
bility phase was followed until Tc = 5.4 x 10® K 
and Tc = 6.9 x 10® K for S300C and S500C re¬ 
spectively, which are good candidates for super- 
collapsar progenitors. Although these end points 
still did not reach the pre-collapse stage (Tc ^ 
10^® K), we did not continue the calculation be¬ 
cause of a numerical difficulty mainly caused by 
very small timesteps encountered at such high 
temperature at the center. For the other model 
sequences, the calculation was terminated imme¬ 
diately before/at the onset of the pair instability 
(Tc = 5 X 10® • • • 10® K), beyond which the evolu¬ 
tionary time is too short to have any further sig¬ 
nificant redistribution of angular momentum. We 
find that the chemically homogeneous evolution 
(CHE) does not occur in any of our model se¬ 
quences (see YDL12 for a detailed discussion on 
the condition of the CHE). 

As shown in Fig. [3l most of the ZAMS mod¬ 
els of the considered mass range have a sufficient 
amount of angular momentum to produce a col¬ 
lapsar, such that any layer of them could form an 
accretion disk around the black hole that would be 
made below it. As mentioned above, these stars 
on the ZAMS rotate as a solid body. As hydrogen 
burning in the core develops, differential rotation 
across the boundary between the contracting con¬ 
vective core and the expanding radiative envelope 
would be created without any transport of angular 
momentum. However, they undergo rapid redistri¬ 
bution of angular momentum mostly via convec¬ 
tion and Eddington-Sweet circulations as long as 
the chemical stratification inside stars is not sig¬ 
nificant. Convection leads to rigid rotation in the 
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Table 1 

Physical properties of the ZAMS models 


Mi 

Tetf 

logL/L© 

F 

'^crit / 

log J 

3 

k 

[M©] 

[10^ K] 




[erg s] 

[10^® cm^ s“b 


300 

106.6 

6.8 

0.774 

0.48 

54.74 

0.93 

0.134 

500 

108.9 

7.1 

0.836 

0.40 

55.05 

1.13 

0.130 

1000 

110.7 

7.5 

0.885 

0.32 

55.46 

1.46 

0.122 

2000 

111.3 

7.8 

0.927 

0.25 

55.86 

1.83 

0.114 

3000 

111.3 

8.0 

0.954 

0.22 

56.08 

2.04 

0.110 

4000 

111.2 

8.1 

0.964 

0.19 

56.24 

2.19 

0.107 

5000 

111.1 

8.2 

0.970 

0.17 

56.36 

2.31 

0.104 

7000 

110.7 

8.4 

0.978 

0.15 

56.54 

2.47 

0.100 

10000 

110.3 

8.5 

0.986 

0.12 

56.70 

2.50 

0.096 

20000 

108.4 

8.9 

0.996 

0.06 

56.88 

1.92 

0.085 


Note. —Each column has the following meaning: Mi: initial mass, Tgff: effective temper¬ 
ature in units of 10^ K, log L/Lq: surface luminosity in units of the solar luminosity, F: 
Eddington factor at the surface, Ucrit/tJK: the critical rotation for the given F at the equa¬ 
torial surface in units of the Keplerian value with rigid rotation, log J: the total angular 
momentum at the critical rotation, j: the specific angular momentum at the critical rota¬ 
tion, k: the dimensionless radius of gyration at the critical rotation (i.e., J — kvcvitMR). 






Fig. 1.— Specific angular momentum distribution in the Pop III models on the ZAMS for 2000, 3000, 
5000, and 10000 M©, as a function of the mass-coordinate (solid line), compared to the specific angular 
momentum at the last stable orbit around the Kerr black hole that would form with the given mass and 
angular momentum (dashed line). 
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Table 2 

Physical properties of the evolutionary models 


Seq. No. 

Mi 

Vi 

j'i/t'c 

log Ji 

Tc.f 

Mf 

Teft.f 

Rf 

log Jt 

Mh© 

GHe 


[Mq] 

[km s-i] 


[erg s[ 

[10®K[ 

[Mo] 

[K] 

[R-q] 

[erg s] 

[Mq] 


S300A 

300.0 

274.0 

0.2 

53.9 

1.0 

298.5 

4380 

5563.7 

53.7 

153.8 

0.6 

S300B 

300.0 

536.8 

0.4 

54.2 

1.0 

291.5 

4081 

5179.2 

53.8 

158.3 

0.9 

S300C 

300.0 

1276.2 

0.9 

54.6 

5.4 

264.0 

4271 

4769.8 

53.8 

200.9 

1.0 

S300M 

300.0 

1276.2 

0.9 

54.6 

0.6 

258.1 

4442 

5642.9 

53.2 

189.7 

0.1 

S500A 

500.0 

454.5 

0.4 

54.4 

0.9 

493.0 

4458 

7530.7 

54.3 

254.4 

0.6 

S500B 

500.0 

598.6 

0.5 

54.6 

0.9 

492.1 

4348 

7912.7 

54.3 

258.6 

0.9 

S500C 

500.0 

1110.8 

0.8 

54.9 

6.9 

464.9 

4157 

7022.7 

54.4 

276.8 

1.0 

S700A 

700.0 

1110.6 

0.8 

55.1 

0.9 

627.9 

4508 

8650.9 

54.6 

417.8 

0.9 

S700B 

700.0 

1384.1 

1.0 

55.2 

0.9 

615.8 

4135 

8246.1 

54.6 

376.9 

0.9 

SIOOOA 

1000.0 

1166.2 

0.9 

55.3 

1.1 

868.0 

3908 

11041.9 

54.4 

492.2 

0.6 

S2000A 

2000.0 

1115.5 

1.0 

55.7 

0.5 

1758.2 

5353 

10291.0 

54.7 

1250.0 

0.3 


Note.—E ach column has the following meaning: M^: initial mass, Vi\ initial rotational velocity at the equatorial surface, 
Vi/vc' ratio of the initial rotational velocity to the critical velocity, log J,: initial total angular momentum, Tc,f: central 
temperature at the end of calculation, M^-. final mass, Teffj: final effective temperature, Rf. final radius, log Jf: final 
total angular momentum, Mhs: mass of the helium core (i.e., the total mass below the hydrogen envelope) at the end of 
the calculation, aHe^ Kerr parameter of the black hole that the entire helium core would make as a result of collapse. 


convective core on the dynamical timescale, and 
Eddington-Sweet circulations in chemically homo¬ 
geneous layers occur on the thermal timescale that 
is much shorter than the evolutionary timescal^. 


For example, in the sequence S300C, solid-body 
rotation is maintained fairly well until the cen¬ 
tral helium mass fraction reaches about 0.75 (The 
top-left panel of Fig. El). The star reaches the crit¬ 
ical rotation soon after the onset of core hydrogen 
burning because of this rapid transfer of angular 
momentum and the gradual increase in the radius 
and surface luminosity, and therefore loses angular 
momentum by the resulting mass loss. However, 
the angular momentum transfer from the core to 
the envelope via Eddington-Sweet circulations is 
significantly slowed down when the chemical strat¬ 
ification becomes strong across the boundary be¬ 
tween the core and the envelope (the so-called 
^—barrier: Mevnet fc Maed^ll997l l. Angular mo¬ 
mentum can be effectively trapped in the core in 
this way and the degree of differential rotation be¬ 
tween the core and the envelope becomes gradu¬ 
ally stronger (Fig. EJ- The final model of S300C 
retains enough angular momentum to produce a 
collapsar (Fig. |3|). 

In the corresponding magnetic model sequence 
S300M, magnetic torques via the TS dynamo are 


^We find that the role of the secular shear instability and 
the GSF instability is minor in our models. 


strong enough to overcome the ^—barrier and 
near-rigid rotation is maintained until the end of 
the main sequence (Top-right panel of Fig. EJ- As 
a natural consequence, this star loses more angu¬ 
lar momentum than in S300C (Bottom panels in 
Fig. [21). In contrary to the corresponding non¬ 
magnetic model, the final model of S300M does 
not retain enough angular momentum for a collap¬ 
sar to occur (see YDL12 for more comprehensive 
discussion on magnetic Pop III star models). 

The initial mass is another key factor to deter¬ 
mine the final outcome. The envelopes of more 
massive stars expand more rapidly both on the 
main se quence and during the post-main sequ ence 
phases ( Marieo et al.ll2003 : lYoon et al.l 2012 ) and 
the amount of angular momentum carried away by 
mass loss, which is proportional to i?^, becomes 
larger. In consequence. Fig. [3] indicates that the 
conditions for collapsar at the final stage can be 
fulfilled only for M, < 700 Mq. From this result, 
we conclude that the initial mass range for poten¬ 
tial super-collapsar progenitors is 300 Mq < M; < 
700 Mq. 


5. Implications for supernovae and gamma- 
ray bursts from super-collapsars 

All of our very massive Pop III star models 
become a red supergiant (RSG; Table 2). It is 
well known that a jet produced in the stellar core 
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-Y =0.93 
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Fig. 2.— Top panels: Evolution of angular velocity on the main sequence in S300C (left panel) and S300M 
(right panel), respectively. The label Yc marks the central mass fraction at the given evolutionary time of 
each profile. Bottom panels: The corresponding integrated angular momentum {J{Mr) = j{m)dm) 
profiles divided by Mr^^. (The evolution of the angular momentum inside stars is better visualized with 
J{Mr)/Mr^^ than with J{Mr) as discussed bv iHeger et al.l ( 2000ll .l The thin contour lines denote levels of 
constant angular momentum in units of erg s. 
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M,/Mo M^Mo 

Fig. 3.— Mean specific angular momentum over the shells. The initial and final distributions of specific 
angular momentum are given by the dotted and solid lines, respectively. The dashed line denotes the specific 
angular momentum for the last stable orbit around a Kerr black hole having all mass and angular momentum 
below the given mass coordinate. 



































via a collapsar cannot easily penetrate a RSG 
envelope and therefore our models may not be 
good progenitors of gamma-ray bursts, unfortu¬ 
nately. To investigate the possibility of a gamma- 
ray burst from our super-collapsar progenitors in 
more detail, we made an estimate on the accre¬ 
tion rate from the disk that the rapidly rotat¬ 
ing layers would make around the black hole, for 
the last models of S500A and S500C as shown in 
Fig- 111 Here we follow ed the approximation by 


Wooslev fc Hegen (120121) : 


M = 


2Mr 

Tff 


— ) 

p- pj 


( 3 ) 


where p = 3 Mr/( 47 rr^) and ts = l/i/247rGp. 
Although the last model of S500C is fairly close 
to the pre-collapse stage, the density of the core 
would further increase until the formation of the 
black hole, and the accretion rate from the core 
material given in the figure should be considered 
a lower limit. The sequence S500A was terminated 
at the onset of carbon burning (see Table 2), which 
is still far from the pre-collapse stage, but the core 
does not retain enough angular momentum for col¬ 
lapsar in this case. The accretion rate from the 
envelope in both S500A and S500C should not be 
affected by this uncertainty, because the hydrogen 
envelope has a much larger dynamical timescale 
than that of the core as shown in the figure. 

The figure indicates that there would be three 
phases of mass accretion via an accretion disk in 
S500C. The first two would be made by accre¬ 
tion from the infalling matters of 3 Mq < < 

47 Mq (assuming the minimum black hole mass 
of 3 Mq) and of 204 Mq < Mr < 277 Mq, and 
the last one from that of the hydrogen envelope 
(379 Mq < Mr < 465 Mq). The expected ac¬ 
cretion rate in the core is very high, implying a 
production of a very powerful jet. However, the 
lifetime of the engine powered by the core mate¬ 
rial would be very short. Most of the core mate¬ 
rial below Mr = 277 Mq would be accreted onto 
the black hole within a second. The outermost 
layer of the core around Mr = 270 — 277 Mq has 
a relatively low accretion rate, but the accretion 
time may not exceed 100 seconds. This is much 
shorter than the jet crossing time through the hy¬ 
drogen envelope (Tcross ^ R*/c ^ 10“^ sec.), posing 
a serious obstacle for the jet penetration and a 
gamma-ray burst (GRB) is not expected from the 


jet produced by the core. 

The most interesting question may be whether 
or not a GRB-like event can be made by the disk 
accretion of the matter in the layers of 434 Mq < 
Mr < 493 Mq and 379 Mq < M^ < 465 Mq of 
the hydrogen envelopes in S500A and S500G, re¬ 
spectively. In particular, the expected accretion 
rate rapidly increases near the surface in Fig. 01 
which res ults from the density i nversion in this 
layer (see IWooslev fc Hegeil 20121 for more discus¬ 
sion on the effect of the density inversion). For 
the case of S500G, a GRB-like event is not likely 
to occur. The energy of the jets from the core in¬ 
jected into the hydrogen envelope would be higher 
than 10^3 erg, if we assume the accretion-to-jet 
conversion efficienc y of about 77 = 10“3 following 
Suwa fc lokal ( 201ll) . Given that Tcross << Tg, the 


forward shock of the jet from the core would be¬ 
come spherical by the time of its breakout from the 
envelope. Because the binding energy of the hy¬ 
drogen envelope is only about 10 ®^ erg, the whole 
hydrogen envelope would be blown away by the 
shock heating, preventing accretion of the outer 
envelope matter onto the black hole. The final 
outcome would be a jet-driven supernova of Type 
HP, rather than a GRB. 


By contrast, for the case of S500A, the inner 
layers are not rapidly rotating, and expected to 
directly collapse to a black hole without making 
a jet. The outermost layers of the hydrogen enve¬ 
lope around the equatorial plane could be accreted 
onto the black hole via an accretion disk to pro¬ 
duce a GRB-like event about two weeks after the 
collapse. In this case, with rj = the gamma- 
ray transient would have L ^ 10^^ erg s“^ lasting 
for about a month. During the last few seconds, 
the accretion rate would dramatically increase as 
shown in Fig. 01 resulting in L ^ 10'^® erg s“^. 

The above discussion is based on the S500A and 
S500G model, but the overall conclusion would 
be the same with the other progenitor models: 
A GRB from very massive Pop HI stars is gen¬ 
erally very difficult to occur because of the red 
supergiant envelope. Gollapsar events in the in¬ 
ner layers would produce a jet-powered type HP 
supernova, if the initial rotation velocity of the 
progenitor were close to the critical limit. In a rel¬ 
atively slowly rotating progenitor, the core would 
not retain enough angular momentum for collap¬ 
sar, but a gamma-ray transient could be produced 
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with the rapidly rotating outermost layers of the 
hydrogen envelope. This event would be marked 
by a long-lasting phase with a relatively low lumi¬ 
nosity (L ^ 10^^ erg s“^ dep ending on ri), which 


may resemble Swif t 1644-1-57 (IBurrows et al.ll2011 
Levan et al.l 201 iL Quataert fc KaseiT 201211 . fol¬ 


lowed by a short phase of a few seconds having 
an enhanced luminosity by factors of 10^ — 10"^ de¬ 
pending on the detailed structure of the layer with 
density inversion, in the rest frame. Such a tran¬ 
sient from redshift of z ^ 20 would be bright in X- 
ray bands and last for a few years in the observer’s 
frame. Given that L > 10®^ erg s“^ is needed for 
a GRB from redshift of abo ut 20 to be detected by 
the X -ray detector BAT (IKomissarov fc Barkov 
2010 ll . the predicted gamma-ray luminosity is too 
low to be observed in near future, unfortunately. 

One uncertain factor in our models is the mass- 
loss rate from very massive Pop 111 stars dur¬ 
ing the red supergiant phase. If the whole hy¬ 
drogen envelope were stripped off before the col¬ 
lapse, favorable conditions for GRB production 
would be more easily fulfilled for the progenitors 
with rapidly rotating cores. As explained above, 
we considered enhancement of mass loss due to 
surface enrichment of heavy elements via chem¬ 
ical mixi ng by extrapolating the mass-loss rate 
given by iNieuwenhniizen fc de Jage^ ( 199fll ). In 
all of our models, the surface metallicity remains 
smaller than about a few times 10“^, and the 
consequent enhancement of the mass loss rate is 
not significant. The mass loss rate might be fur¬ 
ther influenced by pulsational instabilities that 
can ea sily occur in s tars c lose to the Eddington 
limit. Baraffe et al. ( 200lh concluded that mass 
loss induced by pulsation would not be efficient: 
only about 24 Mq is expected to be lost from a 
500 Mq star. This is not significant compared 
to the amount of the lost mass in our models. 
However, the mechanism of RSG winds and its 
dependence on the surface abundance of GNO el¬ 
ements are not well understood yet and remains 
the biggest uncertain factor in our conclusions. 

Another issue that is worth discu ssion is the fi¬ 
nal st ructur e of very massiv e stars . Siiwa fc lokal 
( 2011 1 and Nagakura et al. (2012) showed that 
the GRB jet can breakout a very massive Pop 
III star if it is a blue supergiant (BSG) progen¬ 
itor, for which accretion from the long-lived hy¬ 
drogen envelope to the black hole plays a key 


role (see also IWooslev fc Hegen 120121) . Some pe¬ 
culiar GRBs like GRB 130925A may be well ex¬ 
plained by a progeni t or of this type as recently 
argued bv IPiro et al. ( 2014 ). The relevant ques¬ 
tion is whether or not very massive Pop III stars 
can end its life as a BSG, rather than a RSG. 


The B SG i nodel of 915 Mr:^ used b v ISuwa fc loka 
( 2011 ) and Nagakura et al. ( 2012[) is taken from 


Ohkubo et al.l (|2009l the Y-1 model), who con¬ 


structed this model with rapid mass accretion 
from a 1 Mq protostar. By contrast, in the 
pr esent study a s well as in the previous ones 
bv iBaraffe et al. ( 2001 ). Marigo et al. ( 2003 ) and 
YDL12, all Pop III stars with M > 300 Mq 
become RSGs at the final stage for both non¬ 
rotating and rotating cases, unless they lose most 
of their hydrogen envelopes. This result seems ro¬ 
bust given that each group used a different pre¬ 
scription for overshooting, semi-convection and ro¬ 
tation. To further investigate how the final ra¬ 
dius varies with different assumptions on the over¬ 
shooting and semi-convection, we calculated a cou¬ 
ple of models without overshooting and a small 
semi-convection parameter (asemi = 0.04) with 
300 Mq and 500 Mq and compared them with 
our fiducial models. As presented in Table 3, this 
hardly causes any difference in the final radius for 
300 Mq. a significant reduction in the final ra¬ 
dius with no-overshooting is found with 500 Mq 
but it still does not become a BSG. This strong 
tendency to become a RSG might be related to 
the fact that very massive stars are very close to 
the Eddington limit: according to the equation of 
state, more expansion of a gas is needed to con¬ 
sume a given amount of energy input as the role 
of radiation pressure becomes more important. 

The discrepancy between Ohkubo et al.l ( 20091 
in particular the result with their Y sequences) 
and the others is probably related to the rapid 
mass accretion throughout the whole evolution¬ 
ary stages that was assumed in Ohkubo et ah. 
This assumption is, however, subject to many un¬ 
certainties. For example, it is observed that in 
their M-1 sequence, for which mass accretion was 
stopped for M > 2,21 Mq, the radius rapidly in¬ 
creases to ^ 1000 Rq during the post-main se¬ 
quence phase, unlike the case of the Y sequences 
where rapid mass accretion continues until the end 
of the evolution. Note that the M-1 sequence in 
Ohkubo et al. was constructed to consider the 
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radiation feedback from the mass accreting star 
that can gre atly reduce the mass a ccretion rate as 
discussed bv lMcKee fc TanI ( 2008h . The feedback 
effect must become more complicated with rota¬ 
tion. In particular, mass accretion should result 
in accretion of angular momentum and the mass- 
accreting Pop III proto-stars could easily reach the 
critical rotation when the mass grows beyond a 
certain limit. Given that mass accretion would be 
effectively halted once the stellar surface reaches 
the critical rotation, the Ohkubo et al.’s assump¬ 
tion of continuous mass accretion that gave the 
BSG solution is not likely to be valid with rota¬ 
tion. This issue will be addressed in a separate 
paper (H. Lee and S.-C. Yoon, in preparation). 

Our conclusion is still subject to modifications 
with binary interactions. For example, in a close 
binary system, the mass accreting star may not 
be easily rejuvenated but have a relatively small 
core mass compared to the total mass, if mass ac¬ 
cretion occurs near the end of the main sequence 
or during the post-main sequence phase. This can 
often make the star remain blue during the late 
evolutionary stages (IPodsiadlowski fc Joss 19891: 


Braun fc Langer 1995[| . The evolution of massive 

binary Pop III stars is thus an interesting topic of 
future work. 


6. Conclusions 


We discussed the possibility of super-collapsars 
using the evolutionary models of very massive 
Pop III stars. We find that Pop III stars with 
Mi > 3000 Mq cannot be born with enough an¬ 
gular momentum for collapsars in the core. This 
is mainly because these very massive stars are 
very close to the Eddington limit for which case 
the critical rotation velocity at the equatorial sur¬ 
face is severely limited to a value significantly be¬ 
low the Keplerian velocity. The potential progen¬ 
itors for super-collapsars that can lead to ener¬ 
getic GRBs should have the initial mass range of 
300 M0 < Mi < 3000 Mq. 

We have to consider the evolution of these 
stars to get a more realistic mass range for 
super-collapsar progenitors. High angular mo¬ 
mentum and strong magnetic fields of large 
scales at the pre-collapse stage are the two es- 
sential conditions needed fo r super-collapsars 
( Komissarov fc Barkov! 2010l: Meszaros fc Rees 




M./Mj, 


Fig. 4.— The expected mass accretion rate onto 
the black hole (solid line) given by Eq. ([3]) as a 
function of the mass coordinate in the last mod¬ 
els of S500A (upper panel) and S500C (bottom 
panel). The dashed line denotes the correspond¬ 
ing free-fall timescale {rg = l/v^247rGp, see the 
text). The rapidly rotating layers for which the 
specific angular momentum is higher than the crit¬ 
ical value for the formation of a Keplerian disk (see 
Fig. |3]) are marked by the color shading. 
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Table 3 

Final radius of very massive stars for different mixing parameters 


Mi 

'^crit /’^K 

/over 

Q^semi 

Mf 

Teft.t 

Rf 

[Mq] 




[Mq] 

[K] 

[R-g] 

300 

0.9 

0.000 

0.04 

251.9 

3821 

4761.8 

300 

0.9 

0.335 

1.00 

264.0 

4271 

4769.8 

500 

0.8 

0.000 

0.04 

469.8 

7293 

2550.6 

500 

0.8 

0.335 

1.00 

464.9 

4157 

7022.7 


Note. —Each column has the following meaning: Mi: initial mass, 
^crit/'i^K- the critical rotation for the given T at the equatorial sur¬ 
face in units of the Keplerian value with rigid rotation, /over: over¬ 
shooting parameter. I.e., the convective layer is extended by /over 
times local pressure scale heights beyond the convectively unstable 
region , asemi: semi-convection parameter as defined bv ILanger et all 
lll983f) . Aff: final mass, Teffj: final effective temperature, Rf. final 
radius. 


201Cll) . Given that magnetic torques can lead to 
efficient braking of the stellar core, these two con- 
diti ons may not be ea sily satisfied simultaneously 
(cf. lYoon et alJl20ldl . Our evolutionary models 
indeed show that the angular momentum condi¬ 
tion can be fulfilled only when we ignore mag¬ 
netic torques. If we only consider hydrodynamic 
processes for the transport of angular momen¬ 
tum like Eddington-Sweet circulations, the angu¬ 
lar momentum conditions for super-collapsars at 
the final evolutionary stage can be fulfilled for the 
initial mass range of 300 Mq M\< 700 Mq. 

However, these stars become red supergiants 
at the pre-collapse stage and therefore produc¬ 
tion of an energetic gamma-ray burst from super- 
collapsar events seems difficult. If the core can 
retain enough angular momentum to produce rel¬ 
ativistic jets, it may lead to a jet-powered type 
IIP supernova. If the initial rotational velocity is 
relatively low for the considered mass range, only 
the outermost layers of the hydrogen envelope may 
have specific angular momentum above the critical 
limit for collapsar. The consequent jet formation 
would produce an ultra-long (about a month), rel¬ 
atively faint {L ^ 10"^^ erg s“^) gamma-ray tran¬ 
sient that is marked by a short spike in the gamma- 
ray luminosity (L ~ 10^® — 10®^ erg s“^) during 
the last few seconds, in the rest frame. 

This work was supported by Basic Science Re¬ 
search (2013R1A1A206I842) program through the 
National Research Foundation of Korea (NRF). 
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